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Abstract — We design low-complexity error correction coding 
schemes for channels that introduce different types of errors 
and erasures: on the one hand, the proposed schemes can 
successfully deal with symbol errors and erasures, and, on the 
other hand, they can also successfully handle phased burst errors 
and erasures. 

Index Terms — Decoding, generalized Reed-Solomon code, 
Feng-Tzeng algorithm, phased burst erasure, phased burst error, 
symbol erasure, symbol error. 



I. Introduction 

Many data transmission and storage systems suffer from 
different types of errors at the same time. For example, 
in some data storage systems the state of a memory cell 
might be altered by an alpha particle that hits this memory 
cell. On the other hand, an entire block of memory cells 
might become unreliable because of hardware wear-out. Such 
data transmission and storage systems can be modeled by 
channels that introduce symbol errors and block (i.e., phased 
burst) errors, where block errors encompass several contiguous 
symbols. Moreover, if some side information is available, say 
based on previously observed erroneous behavior of a single 
or of multiple memory cells, this can be modeled as symbol 
erasures and block erasures. 

In this paper, we design novel error correction coding 
schemes that can deal with both symbol and block errors and 
both symbol and block erasures for a setup as in Figure [T] 

• Every small square corresponds to a symbol in F = 
GF((7), where q is an arbitrary prime power (In applica- 
tions, q is typically a small power of 2.) 

• All small squares are arranged in the shape of an m x n 
rectangular array. 

• We say that a symbol error happens if the content of a 
small square is altered. We say that a block error happens 
if one or several small squares in a column of the array 
are altered. 

• Similarly, we say that a symbol erasure happens if the 
content of a small square is erased and we say that a 
block erasure happens if all small squares in a column 
of the array are erased. 

We can correct such errors and erasures by imposing that 
the symbols in such an array constitute a codeword in some 
suitably chosen code C of length ran over F. The two main 
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Fig. 1. Array of size m X n with symbol block errors and erasures. Here, 
m = 8, n = 20, and there are symbol en'ors at positions (2, 2), (4, 2), and 
(6, 14), a symbol erasure at position (7, 3), block errors in columns 5 and 8, 
and a block erasure in column 16. 



ingredients of the code C that is proposed in this paper are, on 
the one hand, a matrix Hin of size m x (mn) over F, and, on 
the other hand, a code C of length n over F. Namely, an array 
forms a codeword in C if and only if every row of the array is 
a codeword in C once the n columns have been transformed 
by n different bijective mappings F™ F™ derived from 
the matrix iJjn. The resulting error-correcting coding scheme 
has the following saUent features: 

« It can be seen as a concatenated coding scheme, however 
with two somewhat non-standard features. First, multiple 
inner codes are used (one for every column encoding), 
and, second, all these inner codes have rate one. 

• One can identify a range of code parameters for C 
for which (to the best of our knowledge) the resulting 
redundancy improves upon the best known. 

• One can devise efficient decoders for combinations of 
symbol and block errors and erasures most relevant in 
practical applications. In particular, these decoders are 
more efficient than a corresponding decoder for a suitably 
chosen generalized Reed-Solomon (GRS) code of length 
mn over F, assuming such a GRS code exists in the first 
place. (Finding efficient decoders for the general case is 
still an open problem.) 

A. Paper Overview 

The paper starts in Section |ll] by considering a simplified 
version of the above error and erasure scenario and of the 
above-mentioned code construction. Namely, in this section 
we consider only block errors and erasures, i.e., no symbol 
errors or erasures. Moreover, an to x ti array forms a codeword 
if and only if every row is a codeword in some code C of 
length n {i.e., there are no bijective mappings applied to the 
columns); in other words, the array code considered is simply 
an TO-level interleaving of C. Our main purpose of Section HIl 
is laying out some of the ideas and tools that will be used in 
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subsequent sections; in particular, it is shown how one can take 
advantage of the rank of the error array in order to increase 
the correction capability of the array code. Nevertheless, the 
discussion in Section III-CI may be of independent interest in 
that it provides a simplified analysis of the decoding error 
probability of interleaved GRS codes when used in certain 
(probabilistic) channel models. 

We then move on to Section [nil which is the heart of the 
paper and which gives all the details of the above-mentioned 
code construction and compares it with other code construc- 
tions. Finally, Section |IV] discusses a variety of decoders for 
the proposed codes. 

B. Related Work 

The idea of exploiting the rank of the error array when 
decoding interleaved codes was presented by Metzner and 
Kapturowski in |22| and by Haslach and Vinck in fT2l . (TSl . 
Therein, the code C is chosen to be a linear [n,k,d] code 
over F, and, clearly, any combination of block errors can be 
corrected as long as their number does not exceed (d — l)/2. 
In II22II and IIT2I . it was further assumed that the set of nonzero 
columns in the (additive) m x n error array E over F is 
linearly independent over F; namely, the rank of E (as a 
matrix over F) equals the number of block errors. It was then 
shown that under this additional assumption, it is possible to 
correct (efficiently) any pattern of up to d — 2 block errors. 
Essentially, the linear independence allows to easily locate 
the nonzero columns in E, and from that point onward, the 
problem reduces to that of erasure decoding. A generalization 
to the case where the nonzero columns in E are not necessarily 
full-rank was discussed in [13 1; we will recall the latter result 
in mode detail in Section III-AI 

The case where the constituent code C is a GRS code 
has been studied in quite a few papers, primarily in the 
context where the contents of each block error is assumed 
to be uniformly drawn from F™. In Q, Bleichenbacher et 
al. identified a threshold, {'m/{m+l)){d—l), on the number 
of block errors, below which the decoding failure probability 
approaches as d goes to infinity and n/q goes to 0. A 
better bound on the decoding error probability was obtained by 
Kurzweil et al. ifTSl and Schmidt et al. Il27l . See also Brown et 
al. Is], Coppersmith and Sudan ||6l, Justesen et al. lfT4l . 
Krachkovsky and Lee ifTTl . and Wachter-Zeh et al. 1291 . 

Turning to the main coding problem studied in this paper — 
namely, handling combinations of symbol errors and block 
errors — a general solution was given by Zinov'ev ll32l and 
Zinov'ev and Zyablov lf33]| . using concatenated codes and 
their generalizations. Specifically, when using an (ordinary) 
concatenated code, the columns of the m x n array are set to 
be codewords of a linear [m, k, d] inner code over F, and each 
of these codewords is the result of an encoding of a coordinate 
of an outer codeword of a second linear [n, K, D] code over 
GF((7'^). It follows from the analysis in ll32l and ll33l that any 
error pattern of up to {} symbol errors and r block errors can 
be correctly decoded, whenever 

2i? + 1 < d(D - 2t) . 



Furthermore, such error patterns can be efficiently decoded, 
provided that the inner and outer codes have efficient bounded- 
distance error-erasure decoders. 

Note that (in the nontrivial case) when 1) > 0, the rate of the 
inner code must be (strictly) smaller than 1. This, in turn, im- 
plies that the overall redundancy of the concatenated code has 
to grow (at least) linearly with ??. A generalized concatenated 
(GC) code allows to circumvent this impediment. We briefly 
describe the approach, roughly following the formulation of 
Blokh and Zyablov H . Given a number r of block errors and 
a number d of symbol errors that need to be corrected, let the 
integer sequences 

1 = do < di < • • • < d„ , 
Do > Di> •••>□„ 

satisfy, for every i = 0,1, . . . ,v, 

2^ + 1 < di{Di ~ 2t) (1) 

(thus. Do > 2{t + -d) + 1 and Dy > 2t + 1). For i = 
1, 2, . . . , w, let H; be an r; x m parity-check matrix of a linear 
[m, ki=m—ri, di] over F. We further assume that these codes 
are (strictly) nested, so that Hi_i forms a proper r^-i x m sub- 
matrix of Hi; the matrix formed by the remaining r^ — ri_i rows 
of Hi will be denoted by 9Hi (we formally define ro = 0, along 
with setting Ho to be an "empty" ro x n matrix). Then an m x n 
array over F is a codeword of the generalized concatenated 
code, if and only if the following two conditions hold: 

(Gl) For i = l,2,...,v, the (partial) syndromes of the 
columns with respect to the partial parity-check matrix 
dHi form a codeword of a code of length n over fi^i"'-! 
with minimum distance Di_i. 

(G2) The columns of the array form a codeword of a code of 
length n over with minimum distance D^,. 

(Note that condition (G2) could be incorporated into condi- 
tion (Gl) by extending the latter to i = u + 1, with Hu_|-i 
taken as an m x m (nonsingular) parity-check matrix of the 
trivial code {0}. Ordinary concatenated codes correspond to 
the case where v = 1 and Do is the "minimum distance" (> n) 
of the trivial code.) It follows from ||32| and ll33l that the above 
array code construction has an efficient decoder that corrects 
any pattern of up to r block errors and up to -d symbol errors. 
See also H], US), El, ED, and the survey Q. 

Recently, Blaum et al. ||2l have proposed new erasure- 
correcting codes for combined block-symbol error patterns. 
The advantage of their scheme is having the smallest possible 
redundancy (equaling the largest total number of symbols that 
can be erased) and an efficient erasure decoding algorithm. 
However, the parameters of their constructions are rather 
strongly limited: first, the array size is typically much smaller 
than q (and, in one embodiment, must in fact be smaller than 
log2 q), and, secondly, verifying whether the construction ac- 
tually works for given parameters becomes intractable, unless 
the number of block erasures or the number of symbol erasures 
is very small. 

In Section IIII-CI we will compare our coding scheme with 
the most relevant of the above-mentioned coding schemes. 
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C. Notation 

This subsection hsts the notation that will be used through- 
out the paper More specialized notation will be introduced 
later on when needed. 

For integers a and b, we denote by (a, b) the set of integers 
{a, a + 1, a + 2, . . ., b— 1}, and (b) will be used as a shorthand 
notation for (0, b). Entries of vectors will be indexed starting 
at 0, and so will be the rows and columns of matrices. For a 
vector u E F" and a subset W C (n), we let {u)w be the 
sub-vector (in I) of u that is indexed by W. The support 
of u will be denoted by supp(M). We extend these definitions 
to any m x n matrix E over F, with {E)w denoting the 
771 X \W\ sub-matrix of E that is formed by the columns that 
are indexed by W. Column j of E will be denoted by Ej, 
and supp(£') will stand for the column support of E, namely, 
the set of indexes j for which Ej ^ 0. The linear subspace 
of F™ that is spanned by the columns of E will be denoted 
by colspan(F). 

The ring of polynomials in the indeterminate x over F will 
be denoted by F\x\, and the ring of bivariate polynomials in 
y and x over F will be denoted by F[j/,a;]Q For a nonzero 
bivariate polynomial ^piy, x) — J2i 'Pi{y)x^ F[y, x], we will 
let deg,j. (p{y, x) stand for the a;-degree of ip{y, x), namely, the 
largest index i for which (pi{y) ^ 0. The y-degree is defined 
in a similar manner The notation F'm „(j/,.T) will stand for 
the set of all bivariate polynomials (p{y,x) £ F[y,x] with 
degj, (p{y, x) < 777 and deg^. (p{y, x) < n. For an element 
^ e F, we denote by Tmiy^S.) the polynomial Y.te(m) Cv'- 

With any m x n matrix E = {eh.j)he(m),je(n) over F, we 
associate the bivariate polynomial 

E{y,x)= e^j/x^' 

he{m).]e{n) 

in Fm.niUi x) (namely, the powers of y index the rows and the 
powers of x index the columns). With each column j of E we 
associate the univariate polynomial Ej{y) = ^hG{m) ^h,jV^', 
thus, E{y,x) = Eje(„) Ej{y)x^ . 

II. Simplified Code Construction 

In this section we consider the simplified scenario men- 
tioned in Section II-AI Namely, we consider only block errors 
and erasures, i.e., no symbol errors or erasures. Moreover, an 
7T7 X 71 array forms a codeword of length mn if and only if every 
row is a codeword in some prescribed code C with parameters 
[n, fc, d] {i.e., there are no bijective mappings applied to the 
columns); equivalently, the array code considered is simply an 
777-level interleaving of C. If C is specified by an {n~k) x n 
parity-check matrix H, then the syndrome matrix S is defined 
to be the 77i x (n—k) matrix S = YH^ , where the m x n 
matrix 

y = F + F 

'We prefer the ordering y,x over x,y because the powers of y and the 
powers of x will be associated with, respectively, the rows and columns of 
m X n matrices like E. 



over F represents the read out (or received) message, where 
the 777 X 77 matrix F over F represents the stored (or trans- 
mitted) codeword, and where the m x n matrix F over F 
represents the alterations that happen to F over time (or during 
transmission). Note that our formalism treats erasures like 
errors, with the side information K C (77) telling us their 
location. 

The subsections of this section are structured as follows. 
In Section III-AI we study the error correction capabilities of 
the interleaved array code, where C is any linear [77, fc,(i] 
code over F. Then, in Section III-BI we present an efficient 
decoder for the special case where C is a GRS code. Finally, in 
Section ITl-CI we present an application of the efficient decoder 
of Section III-BI for the probabilistic decoding of the array 
code under the assumption that the block errors are uniformly 
distributed over F™. 

A. Block Errors and Erasures with Rank Constraints 

We start with Theorem |2] below that generalizes the results 
of II22I and lfT2l to the case where the set of nonzero columns 
of the 777 X 77 error array F are not necessarily linearly inde- 
pendent. Note that this theorem was akeady stated (without 
proof) in the one-page abstract [13J for the error-only case 
{i.e., no block erasures). We include the proof of the theorem 
not just for the sake of completeness, but also because the 
proof technique will be useful in Section |III] we well. 

Toward proving this theorem, the following lemma will be 
helpful. 

Lemma 1 Let C be a linear [77, k^d] code over F and let Z 
be a nonzero m x n matrix over F such that each row in Z 
is a codeword of C. Then 

|supp(Z)| - rank(Z) >d-l. 

Proof: Write J = supp(Z), t = \ J\, and ^ = rank(Z), 
and apply the Singleton bound to the linear \t, fi, >d] code 
over F that is spanned by the rows of {Z)j. ■ 

Theorem 2 Let C be a linear [77, fc, d] code over F and let H 
be an (n—k) x 77 parity-check matrix of C over F. Fix K to 
be a subset of (n) of size r. Given any m x (n—k) (syndrome) 
matrix S over F, there exists at most one m x n matrix E 
over F that has the following properties: 

(i) S = EH^ , and— 

(ii) writing K ~ {n) \ K, the values t = |supp ((F);j^)| and 
jjL = rank((F);^) satisfy 

2t + r<d + ^.i-2. (2) 

Proof: We consider first the case where K is empty. The 
proof is by contradiction. So, assume that F and F are two 
distinct 777 x 77 matrices over F that satisfy conditions (i)-(ii). 
Write 

imax = max{|supp(F)|, |supp(F)|} , 
Mmax = max{rank(F),rank(F)} , 
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and define 

J = supp(£') , 
J = supp(i?) , 
Q = supp(i?) n supp(i?) . 

Consider tlie array Z = E — E. By condition (i) we get that 
every row in Z is a codeword of C. Now, 

|supp(Z)| < \J\ + \J\-\Q\ , 
rank(Z) > 

A'' max \Q\, 

and, so, 

|supp(Z)| - rank(Z) < |J| + |J| 

/^max 
— 2tn:^ax ~ Mmax 

< d- 2 , 

where the last inequahty follows from condition (ii). Hence, 
by Lemma [U we conclude that Z ~ 0, namely, E = E, which 
is a contradiction to the initial assumption. 

Next, we consider the case where r ^ \K\ > 0. First, note 
that condition (ii) implies that r<d + iJ, — 2 — 2t<d— t — 
2; in particular, every subset of r columns in H is linearly 
independent. By applying elementary linear operations to the 
rows of H, we can assume without loss of generality that the 
first r rows of {H)k contain the identity matrix, whereas the 
remaining n—k—r rows in {H)x are all-zero. Let H be the 
(n—k—r) X (n— r) matrix which consists of the last n — k — r 
rows of {H)^: the matrix is a parity-check matrix of the 
hnear [n— r, k] code C over F obtained by puncturing C on the 
positions that are indexed by K. Let S be the m x [n—k—r) 
matrix which consists of the last n — k — r columns of S. We 
have 

S - {E)t^H^ . (3) 

Replacing (i) by ^ and C by C, we have reduced to the case 
where K is empty. The result follows by recalling that the 
minimum distance of C is at least d — r. ■ 
The proof of Theorem|2]in ll22l and lfT2l . which was for the 
special case r — and fi — t— |supp(i?)|, was carried out by 
introducing an efficient algorithm for decoding up to t < — 2 
errors. In that algorithm, Gaussian elimination is performed on 
the columns of S, resulting in an m x (n—k) matrix S'P^, 
for some invertible (n—k) x (n—k) matrix P over F, such 
that the first /i columns in SP~^ form a linearly independent 
set while the last n — k — fi columns in SP~^ are all-zero. (As 
a matter of fact, through this Gaussian elimination, one finds 
the value of fi.) From condition (i) we then get that 

SP'^ = E{PHY . (4) 

Let H' be the [n-k-ii) x n matrix that is formed by the 
last n — k — ji rows of PH. It follows from (|4|l that the 
columns of H' that are indexed by supp(£') must be all- 
zero. Furthermore, all the remaining columns in H' must be 
nonzero, or else we would have /i + 1 < d linearly dependent 
columns in H. It follows that supp(£') is the unique subset 
U C (n) of size fx such that {H')ij is all-zero. Once the 
decoder identifies supp(£'), the entries of E can be found 



by solving linear equations. This decoding algorithm can be 
generalized to handle the case where r = \K\ > {), m the 
spirit of the last part of the proof of Theorem |2] replace E, 
H, and S by {E)j^, H, and S, respectively. 

As pointed out in lfT3l . when K is empty and the difference 
t — fi = |supp(i?)| — rank(i?) is assumed to be equal to some 
nonnegative integer b, then the decoding algorithm of ll22l 
and lfT2] can be generalized into finding a subset U C (n) of 
sizet < ((i+/i)/2-l such that colspan((i/')(7) (C F"'-''-'') 
has dimension b. Letting y be a subset of U of size b such 
that rank((i/')y) = b, the subset W = U\V will then point 
at t — b = jj, columns of E that form a basis of colspan(i?); 
these /i columns, in turn, are flagged as i^i erasures. We will 
then be left with |supp(£') \ W\ — t — fi ~ b nonzero columns 
in E which are yet to be located, but these can be found 
by applying to the received array, row by row, any bounded- 
distance error-erasure decoder D for C. Indeed, since 

25 + /i = 2i - ^ < d - 2 , 

such a decoder can uniquely recover E given the set W 
of erasure locations. The extension to r = \K\ > is 
straightforward. 

Note, however, that even when we are allowed to apply the 
decoder V at no computational cost, we do not know how 
to find the subset U efficiently as b becomes large. In fact, 
if b is large and /i is small, we may instead enumerate over 
the set W which indexes a basis of colspan(i?), then use V 
to reconstruct a candidate for E, and finally verify that we 
indeed have rank(£') = rank((i?)i4.') = ^■ 

B. The GRS Case 

In this section, we present an efficient decoder for finding 
the error matrix E under the conditions of Theorem |2] for the 
special case where C is a generalized Reed-Solomon (GRS) 
code (this decoder will then be used as a subroutine in one 
of the decoders to be presented in Section II Vb . Specifically, 
hereafter in this section, we fix C to be an [n, fc, d=n— fc+1] 
GRS code Cqrs over F with a parity-check matrix 

Hgks = (a; ),g<d-i>jG(n> ' 

where ao, ai, . . . , an~i are distinct nonzero elements of F. 

Let E — ieh.j)he{m),je{n) be an m x n (error) matrix over 
F, and let K and J be disjoint subsets of (n) such that 

J C supp(£;) C KUJ 

(the set J indexes the erroneous columns and K indexes the 
erasure locations). For the matrix E, define the syndrome 
array as the following m x {d—1) matrix S (equivalently, the 
bivariate polynomial S{y,x) E F^.d-iiViX)): 

S = EH^^g 

(compare with condition (i) in Theorem |2). In addition, define 
the error-locator polynomial A(a;) and the erasure-locator 
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polynomial M(a;) by 

A(a;) = 11(1 ~a,x) , 

M(x) = n (1 - "j-^) ' 

respectively. Also, let the modified syndrome array be the 
unique m x (d— 1) array a over F (namely, the unique 
bivariate polynomial a{y, x) in Fm.d-iiu, x)) that satisfies the 
polynomial congruence 



cr{y, x) = S{y, x)M{x) (mod x' 



Finally, the (bivariate) error-evaluator polynomial il.{y,x) 



'Ehe(m) ^hix)y'' is defined by 



jeKUJ j'&{KUJ)\{j} 



Lemma 3 Write t = \ J\, r = \K\, and jj, — rank((iJ)j), and 
suppose that 

2t + r<d + fi-2 

(see ©j. Let X{x) be a polynomial in F[x] and Lo{y,x) = 
Tlihe^m) ^h{x)y'^ be a bivariate polynomial (of y-degree less 
than m) in F[y,x] such that the following conditions are 
satisfied: 

(PI) (j{y, x)X{x) = uj{y, x) (mod x"^"^) , and — 
(P2) degA(.T) < (d+/i-r)/2 - 1 and 
deg^a;(y,x) < {d+fi+r)/2-l . 

Then there is a polynomial u{x) G F[x] such that 

X{x) — A{x)u{x) , 
uj{y,x) = n{y,x)u{x) . 

Proof: First, by the key equation of GRS decoding, it 
is known that (P1)-(P2) are satisfied for A(a;) = A(a;) and 
uj{y,x) = il{y,x) (see, for example, |24, Section 6.3 and 
pp. 207-208]). 

Now, let A(a;) and io{y, x) satisfy (P1)-(P2), fix £ to be any 
index in J, and let [/ be a subset of J of size ^ such that 
i € U and rank((£')[/) = /i (recall that J C supp(£^) and, 
so. Eg 0). Let a = {ah)he{rn) be a row vector in F™ such 
that {aE)ij is nonzero on — and only on — position £. 

Let J be the smallest subset of (n) such that 
J C supp(a^;) C A' U J ; 
note that £ G J and that J Q {i} \J (J \ U) and, so, \ J\ < 



t — fj, + 1. Define 

Hx) = Y[il-a,x) , 

^{x) = ^ — 3" • ^ ah^h{x) 



n (1 - Oij"x) 
3"eu\{e} 



he{m) 



jeKuJ j'e{KuJ)\{]} 
a{x) = ^ ahah{x) , 

u}{x) = ^ ahOJhix) . 

Observing that d-{x) is the modified syndrome polynomial 
that corresponds to the row vector aE, we get from the key 
equation of GRS decoding that 



A{x)a{x) = n(x) (mod x"-') 



Moreover, 



?cd(A(x-),f^(a:)) = l, 



(5) 



(6) 



degA(a;) ^ \J\ 



<t-fi+l < — ^ , (7) 



degi)(a;) < |A'| + |J| <r + t-^ + l< ^^^-^ . (8) 

Multiplying both sides of ^ by X{x) we obtain 

A{x)X{x)a{x) = X{x)n{x) (mod x'^-^) . 

On the other hand, by (PI) we have 

X{x)a-{x) = Lj{x) (mod x'^~^) , 

and so, combining the last two congruences, we get 

A{x)lo{x) = X{x)n{x) (mod x"^-^) . 

Now, from (P2) and ©-([Sll it follows that the degrees of the 
products on both sides of the last congruence are less than 
d—1. Hence, this congruence is actually a polynomial equality: 

A{x)lu{x) = X{x)Cl{x) . 

Thus, from (|6]l we get that X{x) is divisible by A(.t); in 
particular, A(a;) is divisible by 1 — a^x. Ranging now over 
all £ in J, we conclude that A(.t) can be written as A{x)u{x) 
for some polynomial u{x) G F[x]. 

Finally, from (PI) and the key equation we obtain 

uj{y,x) = a{y,x)X{x) 

= a{y, x)A{x)u{x) 

= n{y,x)u{x) (mod x'^'^) , 

from which the equality io{y, x) = u{x)il{y, x) follows again 
by computing degrees: from (P2) we get 

d + Li + r 
deg^ uj{y,x) < 1 < d - 1 
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and 



degx {n{y,x)u{x)) 

= ( deg X{x) - deg A(.t)) + deg^ n{y, : 
< deg A(2;) + r 
d + fi + r 



< 



< d-1 



- 1 



This completes the proof. ■ 
Fix J and K to be disjoint subsets of (n), write i = | J| 
and r ~ \K\, and let be an m x n matrix of rank ^ over 
F such that J C supp(S) C K U J. Define 



(9) 



and 



= ((i?GRs)/f)"' . (10) 

Clearly, the syndrome array S that corresponds to E can be 
decomposed into 

S = S^'^ + , 

and after multiplying the respective bivariate polynomials by 
the erasure-locator polynomial M(a;) — Yij^xi^ ~ ajx), we 
get 

a{y,x) = S{y,x)M{x) 

= S^^\y, x) M{x) + S^^\y, x) M{x) (mod x'^''^) . 

(11) 

Now, recall (from the key equation) that the coefficients of 
x^ ,x^'^^ , . . . ,x'^~'^ in S^^''{y,x)M{x) are all zero, which 
means that the respective coefficients in a{y,x) are equal to 
those in S^'^^y,x)M{x). 

Let S denote the m x (d—l—r) matrix 

5'= io'h.i)he{m},ie{r,d~l} ■ 

It follows from ©-(dU that 

S^{E)j{{AMHGBs).jy , (12) 

where is a (d—l—r) x (d— 1) matrix over F whose first 
row consists of the r + 1 coefficients of M(a;) in decreasing 
order (padded with 2 — r zero entries), and each subsequent 
row is obtained from its predecessor by a shift one position to 
the right (compare (fTZt with (|3]l). Hence, any column indexed 
by .1 G J in the (d—l—r) x t matrix {AmHgrs)j takes the 
form 



[1 a J 



■ ■ a. 



(where M{a~^) =/= since K (1 J ^ 0). Hence, under the 
assumption that t < d — 1 — r (which, in fact, holds whenever 
2t + r < d + - 2) we get that rank((AMi?GRs) j) = t. It 
now follows from (fTZt that 

colspan((i?)j) = colspan(S') (13) 

(provided that t < d — 1 — r); equivalently, for every a G F"\ 

a{E)j = <^ aS = . 



Input: 

• AiTay Y of size m x n over F. 

• Set K of size r of indexes of column erasures. 
Steps: 

1) Compute the m x (d— 1) syndrome array 

5 = YHlns ■ 

2) Compute the modified syndrome array to be the unique m x 
(d— 1) matrix a that satisfies the congruence: 

a{y,x) = S{y,x)M{x) (mod a;''""') , 

where 

Mix) = n (1 - "^■^) ■ 

Let fj, be the rank of the m x (d—l — r) matrix S formed by 
the columns of a that are indexed by (r, d — 1). 

3) Using the Feng-Tzeng algorithm, compute a polynomial X{x) 
of (smallest) degree A < (d+/i— r)/2 such that the following 
congruence is satisfied for some polynomial u}(y,x) with 
deg^a;(y,a;) < r + A: 

a{y,x)\{x) = uj{y,x) (mod a;'*^"') . 

If no such \{x) exists or the computed \(x) does not divide 
nj6(n>(l ~ ctjs;) then declare decoding failure and Stop. 

4) Compute the m x n error array E by 



EAy) = { 



\ ^ ' [ if j £ K 



A(q-^).M'(c«-^) 

, otherwise 

where (•)' denotes formal differentiation. 
Output: 

• Decoded array Y — E of size m x n. 
Fig. 2. Decoding of interleaved GRS codes. (See Section Hl-Bh 



In particular, rank((i?),7) = rank(5). In fact, every m x s 
sub-matrix of S which consists of s > t consecutive columns 
of S has the same rank as {E)j. 

Given cr{y,x) — J2he{m) '^hi-'^)y'^ number of 

erasures r, we can use the Feng-Tzeng algorithm ||9l to 
find efficiently a polynomial X{x) ~ 'J2iLa'^i^^ in 
of smallest degree A such that (PI) is satisfied for some 
w(y,a;) = 'Ehe{m)^hix)y'' where deg^w(y,a;) < r + A. 
In other words, for every h G (m), the sequence (cr/i.OitTr 
satisfies the linear recurrence 



A 

i=0 



= 0, r + A<j<d-2. 



Under the assumption that 2t + r < d + /i — 2, we get from 
Lemma [3] that the polynomial A(a;) found equals the error- 
locator polynomial h.{x) (up to a normalization by a constant). 
From the roots of \{x) one can then recover the set J. 

The decoding algorithm is summarized in Figure |2] Next, 
we analyze its complexity. The syndrome computation step 
(Step [T] in the figure) can be carried out using 0{dmn) 
operations in F. Step |2] requires 0{drm) = 0{d^m) oper- 
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ations to compute a{y, x) and 0{dra ■ min((i, m)) operations 
to compute the rank ji. The application of the Feng-Tzeng 
algorithm in Step [3] requires 0{(Pm) operations, and, finally. 
Step |4] requires 0{dn) operations for the Chien search and 
0{(Pm) operations for computing the nonzero columns of 
E. Overall, the decoding complexity amounts to 0{dmn) 
operations for syndrome computation, 0{dn) for the Chien 
search, and O(d^m) for the remaining steps. 

We note that Lemma |3] and the decoding algorithm in 
Figure |2] apply also to the more general class of alternant 
codes over F, with d now standing for the designed minimum 
distance of the code. Specifically, we apply the lemma and 
the decoding algorithm to the underlying GRS code over the 
appropriate extension field of F: that GRS code has minimum 
distance d and contains the alternant code as a sub-field sub- 
code. 

C. Application to Probabilistic Decoding 

We next provide an application of the efficient decoding 
algorithm of Figure |2] for the following channel model: an 
m X n transmitted array over F is subject to at most a 
prescribed number t of block errors (and no erasures), such 
that the value {i.e., contents) of the block error in each affected 
column is uniformly distributed over F™, independently of the 
other block-error values. (Note that in this formulation of the 
channel, there is a probability of 1/ that an affected column 
will in fact be error- free; we have elected to define the channel 
this way in order to simplify the analysis.) 

We consider the decoding problem given that the m x n 
transmitted array belongs to the code described in Section ITl-BI 
namely, an rn-level interleaving of an [n,k,d=7i—k+l] GRS 
code over F. Let J be the index set of the columns that were 
affected (possibly by an error-free block error), where \J\ < t, 
and let p be the rank of the error array E. The decoder of 
Figure |2] will fail to decode (or will decode incorrectly) only 
when the inequality in Lemma |3] is violated, namely, only 
when /i<2|J| — d+l. Under the assumed statistical model 
on the error arrays, it is easy to see that 

Prob{rank(£') < fi} = Prob {rank((i;),7) < p} 

^ q-{m-i^){\J\-,^) . (l + o(l)) 

(see 1201 p. 699]), where o(l) is an expression that goes to 
as g — > oo. Hence, when to > d — 1, the decoding failure 
probability of the algorithm in Figure |2]is bounded from above, 
up to a multiplicative factor 1 + o(l), by 

^-(m+d-l-2|,/|)(ti-l-| J|) ^ ^-(m+d-l-2t)((i-l-t) 

For TO > d — 1, this bound is (considerably) better than those 
given in |I3] and ifTSl . and is comparable to that in ll27l when 
TO is much larger than d. 

III. Main Code Construction 

We come now to the main code construction of this paper, 
namely the code construction C = (C, iJjn) that was outlined 
in Section ID In particular. Section IIII-AI gives all the details 
of the channel model and the code construction. Section IIII-BI 



presents the error and erasure correction capabilities of the 
code, and Section IIII-CI discusses a variety of embodiments 
based on specific choices for the code C and the matrix 
-ffin, and compares them with alternative code constructions. 
(Decoding algorithms for C will be discussed in Section HVl) 

A. Channel Model and Code Definition 

We consider the following channel model. An m x n array 
r over F is stored (or transmitted), and F is subject to 
the following error and erasure types (see also Table |I] and 
Figure [TJ: 

(Tl) Block errors: a subset of columns in F that are indexed 

hy J C {n) can be erroneous. 
(T2) Block erasures: a subset of columns in F that are indexed 

by C (n) \ J can be erased. 
(T3) Symbol errors: a subset of entries in F that are indexed 

by £ C (m) x {{n)\{K.lJ J)) can be erroneous. 
(T4) Symbol erasures: a subset of entries in F that are indexed 

by 7?. C ((to) X ((n) \ /C)) \ C can be erased. 
Let the m x n matrix £ over F represent the alterations that 
happen to F over time (or during transmission). Then the read 
out (or received) message is given by the m x n matrix 

T = T + £ 

over F. Note that our formalism treats erasures like errors, 
with the side information JC and TZ telling us their location. 
(Of course, the sets J' and C are not known to the decoder a 
priori.) 

Write T = |J|, p = |/C|, t) = \C\, and g = |7^|. The 
total number of symbol errors (resulting from error types (T 1) 
and (T3)) is at most mr + 1? and the total number of symbol 
erasures (resulting from erasure types (T2) and (T4)) is at 
most mp + g; hence, we should be able to correct all error 
and erasure types (T1)-(T4) (when occurring simultaneously) 
while using a code of length mn over F with minimum 
distance at least m(2T + p) + 2i9 + ^) + 1. However, such 
a strategy does not take into account the fact that errors of 
type (T1)-(T2) are aligned across the to rows of the to x n 
array F. The next construction is designed to take advantage 
of such an alignment. 

Definition 4 Let C be a linear [n, k,d] code over F, and let 
i/in be an m x (mn) matrix over F that satisfies the following 
two properties for some positive integer 5: 

(a) Every subset of d — 1 columns in Hi^ is linearly inde- 
pendent (namely, iJjn is a parity-check matrix of a linear 
code over F of length mn and minimum distance at least 
S), and — 

(b) writing 

Hin = ( -ffo I ^^1 I • • ■ I Hn-i ) , 

with Hq, Hi, . . . , Hn-i being m x m sub-matrices of 

Hin, each Hj is invertible over F. 
Given C and Hi^, we define C = {C,Hin) to be the linear 
[mn, mk] code over F which consists of all m x n matrices 

F = ( Fo I Fi I ... I r„_i ) 
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TABLE I 

Types of errors and erasures under consideration. 





error 


erasure 


block 


(Tl) 

column set J 

\J\=r 


(T2) 

column set K 

\K\ = p 


symbol 


(T3) 

location set £ 

\C\ = •& 


(T4) 

location set TZ 

\n\ = g 



over F (where Tj stands for column j ofT) such that each 
row in 



z = ( iJoro 

is a codeword of C. 



HiFi I ... I ff„_ir„_i ) (14) 

□ 



One can view the code C as a (generalized) concatenated 
code, where the outer code is an m-level interleaving of C, 
such that an m X n matrix 



Z=( Zo\ 



Zn-l 



over F is an outer codeword if and only if each row in Z 
belongs to C. Each column in Z then undergoes encoding by 
an inner encoder of rate one, where the encoder of column j 
is given by the bijective mapping Zj i-^ H~^Zj. 

B. Error Correction Capabilities 

This subsection discusses what combinations of errors and 
erasures of the types (T1)-(T4) can be handled by the code 
C = (C,i?in) that was specified in Definition ID 

Theorem 5 There exists a decoder for the code C that 
correctly recovers the transmitted array in the presence of 
errors of types (T1)-(T4) (which may occur simultaneously), 
whenever r \ J\), p |/C|), d \C\), and g (= |7^|) 
satisfy 

2t + p<d-2 , 
2i} + g<S-l . 

Proof: Using puncturing as in the proof of Theorem |2] it 
suffices to prove the theorem for the case where K, is empty, 
i.e., there are no erasure events of type (T2). 

The proof is by contradiction. So, assume that £ and £ are 
two distinct mx n matrices that correspond to error events of 
types (Tl), (T3), and (T4), with the respective sets J', J , C, 
and C satisfying 

2rn,ax < - 2 , 

Tmax = maxfUl, \JW , 

^?max = max{|£|, |£|} , 



where 



and g = \Ti\ (the symbol erasure set TZ is the same for both 
£ and £). We will have reached a contradiction once we have 
shown that £ — £ is a codeword of C if only if £ = £. 
Let Z be the m x n matrix 

Z = {Ho{£o—£q) I Hi{£i—£i) | . . . | Hn-l{£n-l—£n-l) ) • 

Observe that since the matrices Hq, Hi, . . . , H^^i are all 
invertible over F, we get that a column in Z is zero if and 
only if it is zero in £ — £; in particular, Z = if and only if 
£^£. Write 

Q = 5upp{£ -£)\{JUJ) , 

namely, the set Q indexes the columns of £ — £ that contain 
errors of type (T3) and (possibly part of) the erasures of 
type (T4). For each j E Q, denote by wj the number of 
nonzero entries in £j — £j (that is, wj = |supp((£'j — £j) )|). 
The total number of nonzero entries in {£ — £)q satisfies: 

<\CU CUTZ\ < 2i?,„ax + g<S -1 . 

Consider the respective columns in Z: 

Z,^H,{£,-£,), jeQ. 

Each Zj is a nontrivial linear combination of wj columns of 
-ffiii, and, obviously, for distinct indexes j these combinations 
involve disjoint sets of columns of Hin- Recalling that every 
subset of X^jeg'^^j (— i5 ~ 1) columns in iJjn is linearly 
independent, we thus get that the set of columns of {Z)q is 
linearly independent, that is, 

rank(Z) > rank((Z)Q) = \Q\ . 

On the other hand, 

|supp(Z)| < \J\ + \J\ + \Q\ < 2w + |Q| < d - 2 + |Q| 
and, so, 

|supp(Z)| - rank(Z) <d-2 . 

Hence, we conclude from Lemma[T]that each row in Z belongs 
to C only if Z ~ 0. Equivalently, £ — £ belongs to C if only 
if £ = £, as promised. ■ 

Remark 6 We draw the attention of the reader to the con- 
dition on r and p in Theorem |5] namely, that the expression 
2r + p be at most d — 2, rather than the (more common) 
requirement that it be at most d—1. It is this slightly stronger 
condition that, implicitly, provides the required redundancy for 
correcting the (additional) symbol errors and erasures. □ 

Remark 7 Theorem |5] indirectly implies a dependence of 
the correction capability of errors of type (T3)-(T4) on the 
parameter to, which, in turn, is part of the specification of 
the errors of type (T3)-(T4). Specifically, the largest possible 
value for S can be the minimum distance of any linear code 
of length mn and redundancy m over F. Of course, one may 
re-arrange the array by grouping together non-overlapping sets 
of s columns, for some integer s > 0, to form an to' x n' array 
where to' = sm and n' ^ n/s (assuming the latter ratio is an 
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integer). The block errors and the symbol errors wiU remain 
so also in this modified setting, except that the block errors 
will be more structured than just being phased with respect 
to the (new) parameter m'. If this additional structure is not 
taken into account, the code C is bound to be sub-optimal. □ 

Remark 8 In contrast to the previous remark, if m is (much) 
larger than the redundancy needed from a linear code of 
length mn and minimum distance S as in Theorem |5] we 
may partition each column in the original array into s new 
columns, thereby forming an m' x n' array, where m' ^ m/s 
and n' = sn, such that m' is (just) the redundancy required 
from a linear code of length mn and minimum distance S. 
Of course, this will increase the number of block errors by a 
factor of s, yet as long as n' is sufficiently small to allow us to 
use a maximum-distance separable (MDS) code for C (as will 
be the case in the examples in Section IIII-CI ). we will incur a 
gain in the redundancy: it will reduce from (2r + p)m + m 
to (2r + p)m + m'. Thus, the code C is suitable for (the 
rather practical) scenarios where the number of symbol errors 
is relatively small compared to the block-error length m. □ 

C. Examples 

In this subsection, we consider various special choices for C 
and Hin and demonstrate the properties of the resulting code C 
(Examples l9l through [T2]|. We also demonstrate the advantages 
of the code C over existing alternatives for handling errors of 
types (T1)-(T4) (Examples [T3] through [T6]l. 

Example 9 We consider first the special case where mn < 
q + 1. Here, we can take C to be an MDS code over F and Hj^ 
to be a parity-check matrix of an MDS code over F. Under 
such circumstances we have d = ?i — fc + 1 and S = m + 1, 
which means that it suffices that the sizes r, p, d, and g satisfy 

2T + p<n-k-l, 
2^ + g<m . 

The redundancy of C, being m{n — k), is then the smallest 
possible for this correction capability: since the total number 
of symbol errors is mr + d and the total number of symbol 
erasures is mp+ g, by the Reiger bound ( |fT9l , ||231 ) we need 
a redundancy of at least m(2r + p) + 21) + g symbols over 
F in order to be able to correct all error types (T1)-(T4). 
Admittedly, the same performance of correction capability 
versus redundancy can be achieved also by a single linear 
[mn, mk] MDS code C over F (which exists under the 
assumption that 77in < q+l). However, as pointed out earlier, 
the use of such a code C does not take into account the 
alignment of error types (Tl) and (T2) across the rows of 
the received mx n array. It is this alignment that allows C to 
achieve the same correction capability using a code C, which 
is m times shorter than C. While we still need for C the 
parity-check matrix Hi^ of an MDS code of length mn, the 
redundancy of the latter code needs to be only m, rather than 
m{n — k). 

To demonstrate the savings that C may offer compared 
to C, consider the simple problem of verifying whether a 



given m X n array T belongs to the code (namely, detecting 
whether errors have occurred). When using C, we will regard 
r as a vector of length mn over F and the checking will be 
carried out through multiplication by an {m{n~k)) x [mn) 
(systematic) parity-check matrix of C, thereby requiring up to 
2m^k{n — k) operations (namely, additions and multiplica- 
tions) in F. In contrast, when using C, we will first compute 
the array Z as in (fT4] | while requiring less than 2m^(n — 1) 
operations in F (one of the matrices Hj can be assumed to 
be the identity matrix); then we will compute the syndrome 
of each row of Z, for which we will need up to 2mk{n — k) 
operations in F. □ 

Example 10 Suppose that mn < q — \ and select C to be an 
[n, k, fc+1] GRS code Cgrs over F as in Section HTB] 

For every j e (n), let 

such that the elements j are distinct and nonzero in F for all 

n G (m) and j e (n); the respective matrix H\n = {Hj)j^i^n) 
is then a parity-check matrix of an [mn, m(7i— 1), m+l] GRS 
code over F. Given any mx n matrix F = (T K,j)Ke{m),je{n) 
over F, the entries of HjTj are given by 

KG (rn) 

Hence, F is in C = (C, iJin) if and only if 

E E rK.,«}/3K.j = , e (m) , ie{d-i). 

(15) 

(Note that if /3k j depended only on k, then C could be seen 
as a two-dimensional shortening of a two-dimensional cyclic 
code; see, for example |26].) □ 

Example 11 We show that sometimes the construction C in 
Example [Tol is an MDS code. Assume therein that m divides 
q — I and that for every j e (n), the multiplicative order of 
aj divides {q — l)/m (thus, each aj has m distinct mth roots 
in F). For every j e (n), select /3o.j,/3i.j, ■ ■ ■ ,/3m-i.j to be 
the distinct roots of aj in F. Then F e C if and only if 

E E rK,,/3:'+"" = , he (m) , leid-i). 

The latter condition, in turn, is equivalent to F being a 
codeword of a GRS code of length mn and redundancy 
{d — l)m over F. □ 

In contrast, the following example shows that sometimes 
C = (C, i/in) is not an MDS code, even when C is MDS and 
is a parity-check matrix of an MDS code. 

Example 12 Suppose that g is a power of 4 and take m = 3 
and n = (q + 2)/3. Select Hi^ to be a parity-check matrix of 
a [q+2, 4] triply-extended GRS code over IMl P- 326] 
and C to be any linear [n, k] code over F. Thus, C is a linear 
[q+2, 3k] code over F. It follows from the already-proved 
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range of the MDS conjecture that C, being longer than q+l, 
cannot be MDS when, say, 2<k<l + ^y/q^. □ 

In fact. Example [12] shows that there are choices of F, n, 
in, T, and for which the construction C = (C,7?in) can 
be realized to correct any r block errors and any -d symbol 
errors, while, on the other hand, there are no codes over F of 
the same length and size as C that can correct any 2Tm + -d 
symbol errors. Examples [13] and [14] below present a larger 
range of parameters where this may happen. 

In Examples [T3] - [T6] we make a running assumption that 
n < q, in which case C can be taken as an MDS code, such as 
a (possibly extended) GRS code. For the sake of simplicity, 
we will consider in these examples only the block-symbol 
error-only case, i.e., no erasures are present. 

Example 13 Given positive r, -d, n, and F = GF{q) (such 
that 2r + 2 < 71 < g), we take Hm to be a parity-check matrix 
of a (possibly extended) shortened BCH code of length mn 
over F, where m is determined by n, and q to satisfy the 
equality 



1 



(2,9-1) • log„(mn) 



(so mn may be larger than q; see |[24l p. 260]). The code C 
is taken as a (possibly extended) [n,k,d] GRS code over F 
where d = 2t + 2. The overall redundancy of C = (C, iJin) 
is then 



(2t + l)m = 2Tm + 1 



q 
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log Amn) 



(16) 



The first term, 2Tm, on the right-hand side of ([T6] l is the small- 
est redundancy possible if one is to correct any r block errors 
of length m. The remaining term therein is the redundancy (or 
an upper bound thereof) of a BCH code that corrects any -d 
symbol errors over F. In comparison, a shortened BCH code 
of length mn over F that corrects any rm + 1? symbol errors 
may have redundancy as large as 



1 



^ • (2(Tm + i9) - 1) 
q 



log (mn) 



It can be verified that the latter expression is larger than (ITSI l 
when mn > g > 4. □ 

Example 14 We repeat the previous example, except that 
we replace the shortened BCH codes with codes that attain 
the Gilbert- Varshamov bound. Denote by reAq{N,D) the 
redundancy of such a linear code of length N and minimum 
distance D over F. Then m is determined by -d, n, and q 
through 

m = redq{mn, 2^? + 1) . 

Taking again C to be an [n, k, d] GRS code over F, the code 
{C,Hia) has redundancy 



while a linear code that attains the Gilbert- Varshamov bound 
and corrects any rm + symbol errors has redundancy 

redq(mn, 2{Tm + d) + 1) . 

Thus, we are interested in the range in which the following 
inequality holds: 

redg(mri, 2{Tm + ■§) + 1) > 2Tm + redq{mn, 2^ + 1) . 

(17) 

As shown in Appendix [A] the inequality dTTI i indeed holds 
whenever 

9-1 



2(Tm + ^ — 1) < mn ■ 



2q - 1 



(18) 



which, in turn, holds whenever 



2t + 1 < n • 



g-1 
2q- 1 



(recall that 2?? < m). 



□ 



Example 15 We compare the performance of C with that of a 
concatenated code C constructed from a linear [m, k, d] inner 
code over F and a linear [n, K, D] outer code over GF(g'^) 
(where n < q). By the Singleton bound, we can bound the 
redundancy of C from below by 

mn- kK > (D-l)m + (d-l)n - (D-l)(d-l) 

= (D-l)(77i + l) -n + d(n-D + l) . (19) 

As we mentioned already in Section II-BI any error pattern of 
up to T block errors and up to -d symbol errors can be correctly 
decoded, whenever 



279 + 1 < d(D-2r) 



(20) 



Since in our setting the values r and -d are prescribed, we can 
minimize (l% over d and D, subject to the inequality ( [20l i. 
Specifically, we define A = D — 2t and, from ( i20l i. we can 
express d in terms of A as d = \ {213 + 1) / A] , in which 
case ([T9] l becomes 



A(777 + 1) + 



A 



2r + 1 



+ (2T-l)(7r7+l)-77 

(279 + l)(?i-2r+l) 



> A(777+ 1) 



A 

2r+ 1) 



- (27? + 1) - (77 - 2r + 1) + (2r - 1)777 
The minimum of ( i22] i over (a real) A is attained at 



(21) 



(22) 



(279 + l)(77-2r + l) 

7T7 + 1 

Yet we need to take into account that both d and A are positive 
integers. 
Case 1: n 

we have Amin < 1 and, so, we take A = 1; namely, for this 
range, it is best to let the inner code handle (exclusively) the 
symbol errors, and then the outer code is left to correct the 
block errors. For this range, the expression (ED becomes 



2t < In this range (of very small 77) 



(2t + 1)777 = 2t777 + red5(77777, 2l? + 1) , 



2r77i + 2'd{n - 2t) , 
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which is never larger than the smallest possible redundancy, 
(2t + l)m, of C. Notice, however, that when H-^^ in C can be 
taken as a parity-check matrix of an MDS code, then to = 21? 
and therefore this range is empty. 

Case 2: n — 2t > (to + 1) • (21? + 1). In this range we 
have Ainin > 219 + 1 and, so, the expression [•] (for d) in (ISTT i 
equals 1; namely, for this range, it is best to regard the symbol 
errors as block errors. Therefore, we take A = 2-d + l and the 
expression jTli becomes 

2(t + '!9)to , 

which is larger than the redundancy of C whenever i9 > 
(assuming that C in C is an MDS code). 

Case 3: < n - 2r < (to + 1) • (2z? + 1). In this range, 
we plug the expression for Amin into (|22] |. resulting in the 
following lower bound expression on the redundancy of C; 

2\/(2i? + l)(n - 2t + l)(m + 1) - (2?9 + 1) 
- (n - 2t+ 1) + (2r - 1)to . 

This lower bound is greater than the redundancy of C when- 
ever 

4(2i? + l)(n-2r + l)(TO+l) 

> (2TO+(27?+l) + (n-2r+l))\ (23) 

Since the left-hand side of (l23l l is a cubic expression while 
the right-hand side is only quadratic, the inequality (l23l l is 
expected to hold for most parameters of interest. □ 

Example 16 In many cases, the redundancy of C is smaller 
even than that of the generalized concatenated (GC) code 
construction defined through conditions (G1)-(G2) in Sec- 
tion II-BI Referring to the notation therein, we first note that 
if Du > 2r + 1, then the contribution of condition (G2) alone 
to the redundancy is already at least (D„ — 1)7ti > (2t + 1)to. 
Hence, we assume that = 2r + 1, in which case, from ([T]i, 
we get that ry > d„ — 1 > 2i?. Thus, condition (G2) induces 
a redundancy of at least 2tto, and condition (Gl) adds a 
redundancy of at least 

V 

^(r, -r,_i)(D,_i-l-2r). (24) 

In Appendix |Bj we show that the expression ( l24l i is bounded 
from below by 

(2i? + 1) In ?? + 27 • + 0(1) , (25) 

where 7 is Euler's constant (approximately 0.5772) ifTOl 
p. 264]. Taking C in C as an MDS code, the redundancy of C 
is then smaller than that of GC codes (with the same correction 
capabilities) whenever i25[ exceeds to. □ 

One may speculate whether Theorem|5]holds for the follow- 
ing more general definition of C; instead of requiring that each 
row of Z in (fT4l i be a codeword of C, require that Z belong 
to a linear [n, k, d] code over GF(g™), where each column of 
Z is regarded as an element of the latter field with respect to 
some basis of that field over F. It turns out that Theorem |5] 



does not hold for this more general setting, as shown by the 
counterexample that is discussed in the following example. 

Example 17 Let d = 2 (i.e., t = p = 0). Assume that n < 
— 1 and let Co, Ci, C2, . . . , C„-i be any n distinct powers 
of an TO X TO companion matrix of some irreducible polynomial 
of degree to over F ll20l p. 106], where Cq = I (the m x m 
identity matrix). Then 

|z = ( Zo I Zi I . . . I Z„_i ) : ^1^1 = o| 

is a linear [n, n—1, 2] code over GF(q™) (with the columns of 
each Z being the codeword coordinates). The respective code 
C would then be written as 

|r = ( To I Ti I . . . I r„_i ) : J2 Cj-^j-r. = o| • 

If (5 > 3 then the 2to columns of Hq and Hi are all nonzero 
and distinct. Therefore, we can select Ci 7^ / so that the 
first column (say) in C\}1\ equals the first column (say) in 
= CqHq. Yet this means that there exists a nonzero F G C 
that contains only two nonzero entries. Therefore, C cannot 
have a decoder that corrects all one-symbol error patterns. □ 

IV. Decoding Algorithms for the 
Main Code Construction 

We now discuss a variety of decoders for the code C = 
(C, Hin) that was specified in DefinitionlH for the case where 
the constituent code C is a GRS code. Section ITV-Al presents a 
polynomial-time decoding algorithm for error and erasures of 
types (T1)-(T4) as far as they are correctable as guaranteed 
by Theorem |5] and provided that the code parameters n, d, 5, 
p, and Q satisfy a certain inequality (see Theorem [T9] below). 
Then, we consider the construction C as in Example [TOl( where 
C is a GRS code and is a parity-check matrix of a GRS 
code) and present some more specialized decoders for this 
construction. Namely, Section IIV-BI discusses a decoder that 
handles errors and erasures of types (Tl), (T2), (T4), but 
not of type (T3), and Section IIV-CI introduces a decoder that 
handles errors and erasures of types (Tl), (T2), (T4), and some 
combinations of errors of type (T3). As of yet, we do not have 
an efficient decoder that corrects all error patterns that satisfy 
the conditions of Theorem |5] 

Assume that to, r, p, d, and q scale linearly with n. If C 
is replaced by a GRS code (if such a code exists) then the 
decoding complexity scales linearly with [n^Y ~ 
of the main purposes of defining the code C = (C, Hi^) and 
formulating the decoders in Sections IIV-BI and IIV-CI is that 
their decoding complexity scales linearly with n^. 

A. Polynomial-Time Decoding Algorithm 

Example [TT] demonstrates one particular instance of the 
construction C = (C, Hin) for which the decoder guaranteed 
by Theorem|5]has an efficient implementation (simply because 
in this case the code C is a GRS code). In this section, we 
exhibit a much wider range of instances for which decoding 
can be carried out in polynomial-time complexity. 
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Specifically, we consider the case where the code C is a GRS 
code over F (and, so, n < g — 1), and Hin is an arbitrary 
m X (mn) matrix over F that satisfies the two properties (a)- 
(b) in Definition |4l The columns of m x n arrays will be 
regarded as elements of the extension field GV{q"^) (according 
to some basis of GV{q™') over F). When doing so, the matrix 
Z in (fT4l i can be seen as a codeword of a GRS code C over 
GF((7™), where C has the same code locators {oij)jtz(n) 
C. 

Let r e C be the transmitted m x ?? array and let T be 
the received m x n array, possibly corrupted by t errors of 
type (Tl) and d errors of type (T3), where r < {d/2) — 1 and 
1? < (5 — l)/2. We first compute an m x n array 



Y = ( i/oTo I i?iTi 



where Y contains t + d < {d + 5 — i)/2 erroneous columns. 
Regarding now y as a corrupted version of a codeword of 
C, we can apply a list decoder for C to Y . Such a decoder 
returns a list of up to a prescribed number L of codewords of 
C, and the returned list is guaranteed to contain the correct 
codeword, provided that the number of erroneous columns 
in Y does not exceed the decoding radius of C. In our 
decoding, we will use the polynomial-time list decoder due 
to Guruswami and Sudan [fTTI (for variation of the algorithm 
that reduce its complexity, see |fT6l and ||30l ). For any GRS 
code of length n and minimum distance d (over any field) 
and for any prescribed list size L, their decoder will return 
the correct codeword as long as the number of errors does not 
exceed \n0L{d/n)'\ — 1, where 0L{d/n) is the maximum 
over s € {1, 2, . . . , L} of the following expression: 



OL.sid/n) 



1 



1 



1 



2{L + 1) 2s V n 
(see lEl Chapter 9.5]). Thus, if L is such that 

neLid/n)>{d + 5-l)/2, (26) 

then the returned list will contain the correct codeword 



Z 



HoTq HiTi 



of C. For each array Z' in the list we can compute the 
respective array in C, 



r' 



-^ri-l^n-l 



and it follows from the proof Theorem |5] that only one such 
computed array F' — namely, the transmitted array F — can 
correspond to an error pattern of up to {d/2) — 1 block errors 
and up to {5 — l)/2 symbol errors. Finding that array can 
be done simply by checking each computed Z' against the 
received array T. 

Remark 18 While the decoding scheme that we have just 
outlined makes essential use of an efficient list decoder for 
the m-level interleaving of C, nothing is assumed about 
beyond properties (a)-(b) in Definition |4] In particular, nothing 
is assumed about the decoding complexity of the code over F 
that is defined by the parity-check matrix H^a. □ 



Our decoding scheme can be generalized to handle also 
erasure of types (T2) and (T4) by applying a list decoder for 
the GRS code obtained by puncturing C on the columns that 
are affected by erasures: this translates into replacing d by 
d — p — g (assuming that the latter value is positive). 

The next theorem characterizes a range of parameters for 
which (|26] | holds for some polynomially-large list size L (and, 
thus, C can be decoded in polynomial time). 

Theorem 19 For C = (C,i7in) such that C is a GRS code 
over F, the decoder guaranteed by Theorem \5\ can be imple- 
mented by a polynomial- time algorithm, whenever 



d- p- g> 2^S{n - p - g) - 5 (27) 

(or, simply, whenever d > 2\fEn — 5 in case p = g = Q). 

Proof: We will assume in the proof that p = g = the 
general case follows by observing that any puncturing of C 
on p + g positions results in a GRS code of length n — p — g 
and minimum distance d — p~ g. 

Our proof will be complete once we identify a poly- 
nomially-large L for which ( |26] | holds. We take L to be such 
that 



> 2 



1 



1 



(28) 



n V ?i n L + 1 

It readily follows from (|27] | that there exists such an L which 
is at most quadratic in n. 
Define s to be 



L- 



L ■ 



1 



From 



we have 

d 

- > — 

n L 



5-1 L+l-s 



s n 

which can also be rewritten as 

'5-1 L 
~2^ ^ 2^ 



d 2s 

- > 

n L — s 



L + 1 



s + 1 
2(1 + 1) 



1 



Multiplying both sides by (L — s)/{2s) and rearranging terms 
yields 

s + 1 L / d\ d 5-^ 

1 r-^ r 1 - - > — 

" " ' 2s V nj ~ 2n 



2{L + 1, 
which is equivalent to 

OLAd/n) > 
This, in turn, implies 



2n 



d + 5-1 
2n 



B. Decoding of Errors and Erasures of Type (Tl), (T2), (T4), 
but not of Type (T3) 

In this section, we present an efficient decoder for the code 
C when constructed as in Example [TO] for the special case 
where i) = \C\ — (no errors of type (T3)). 

An m X n matrix F is transmitted and an m x n matrix 

T = T + £ 



13 



is received, where 8 ~ {£K,j)Ke{ni).je{n) is an m x n error 
matrix, with J' (C {n)) (respectively, K, (C (n))) indexing the 
columns in which block errors (respectively, block erasures) 
have occurred, and TZ (C (m) x (ri)) is a nonempty set of 
positions where symbol erasures have occurred. We assume 
that d, T (= I J'l), and p (= |/C|) satisfy 



2t + p<d-2 
and that g (= \TZ\) satisfies 

< g < m . 

Define 

Y = ( i^oTo I i/iTi I 

and 

^ = i(^h,j)he{m)J£{n) 



(29) 



Clearly, 



Y = Z + E 



where Z is given by ( fT4b . In particular, every row in Z is a 
codeword of Cgrs- 

Next, write TZ — {iKe,je)}ee{g}- Fix some i G (g) and 
define the following g univariate polynomial (of degree g—1) 



it) 



y 



n 



1 - f3.jy 



(30) 

, e^{g) , (31) 



and let e*^^) = (ej^-')jg^„) denote row £i— lofthe {Tn+g—l)xn 
matrix B^^\y)E{y, x). We have 

supp(e(^))C j7U/CU{ja, fG(^>). 

Indeed, the contribution of a symbol erasure at position (k, 
in £ to the column Ej{y) of E{y,x) is an additive term of 
the form 



1 - iP^.jyr 
1 " P^,,y 



(where we recall the definition of Tm(-; •) from Section IFCI i; 
so, if 7^ then the product 



B(')(y)-e„,,- 



B(^)(;/) 



have 



is a polynomial in which the powers y 
zero coefficients. 

Now, every row in the (m+g—l) x n array Z'^^\y,x) = 
'Q^^\y)Z{y, x) is a codeword of Cqrs- Therefore, by applying 
a decoder for Cqrs to row p — 1 of Z^^") with p + 1 erasures 
indexed by /C U {j^}, we should be able to decode the vector 
e'^^\ based on our assumption 



It follows from the definition of e'^^ that for every j G (n). 



(e-i) 



In particular, 



B, 



(1) 



(1) 



B 



\Bo 



(e-i) 



B 



(e-i) 



B 



0-1 



(e-i) 
e-i y 



(32) 



1-1 



Because £ E (g) was arbitrary, we are able to recover the error 
values in £ at the positions TZ. Namely, 

This, in turn, allows us to eliminate the symbol erasures from 
E. 

Figure [3] summarizes the decoding algorithm of a com- 
bination of errors of type (Tl), (T2), and (T4). The com- 
plexity of Step [T] is 0[{d + ?T!,)mrt) operations in F (see 
the discussion that precedes Example [TO] i. Step |2] requires 
0{dpg) operations. Each iteration in Step |3] requires 0{dg) 
operations (for Step|3a), 0{d^) operations (for Step |3b]i. and 
0{dm) operations (for Step |3c), totaling to 0(d{d + m)g) 
for Step|3] Step|4]requires 0{d'^m) operations to compute the 
error-locator and error-evaluator polynomials, and additional 
0{dn) for the Chien search. Finally, Step |5] requires 0{drn?) 
operations. To summarize, the decoding complexity amounts 
to 0(((i + m)mn) operations for syndrome computation, 
0{dn) for the Chien search, and 0{d{d + m)nij for the 
remaining steps. 



C. Decoding of Errors and Erasures of Type (Tl), (T2), (T4), 
and with Restrictions on Errors of Type (T3) 

In this section, we consider the decoding of C when 
constructed as in Example [TO] under certain assumptions on 
the set C, namely, under some restrictions on the symbol error 
positions (errors of type (T3)). These restrictions always hold 
when |£| < 3 and d is sufficiently large. 

Specifically, we consider the case where each column, 
except possibly for one column, contains at most one symbol 
error The general strategy will be to locate the positions 
of these errors, thereby reducing to the case considered in 
Section IIV-BI We use the same notation as in that section, 
except that the set C is not necessarily empty and that (for 
reasons of simplicity) the set TZ is empty. As in Section [IV-BI 
the number r of block errors and the number p of block 
erasures satisfy 2t + p < d — 2. 

When I? = |£| > 0, we write C = and 
assume that that there exists a w E {-d) such that the values 



14 



Input: 

• Array T of size m x n over F. 

• Set IC of indexes of column erasures. 

• Set TZ = {{i^e, je)}e£{e) of positions of symbol erasures. 
Steps: 

1) Compute the m x {d—1) syndrome array 

S = [ HaTo I H\T\ | • ■ • | Hn-iT n-i ) Hgrs ■ 

2) Compute the modified syndrome array to be the unique g x 
(d— 1) matrix a that satisfies the congruence 

a{y,x) = S{y,x)l[{l - a,x) (mod {x^-^y"}) . 
jeic 

3) For every £ £ (g) do: 

a) Compute row g — 1 in the unique g x (d— 1) matrix cr'^' 
that satisfies the congruence 

a^'\y,x)=B'^'\y)aiy,x){l-a,,x) 

(mod {x'^''^,y^}) , 

where B('^'(j/) is as in i30l . 

b) Decode e^^' (i.e., entry in e'*') by applying a decoder 

for Cgrs using row ^3 — 1 in o"'*' as syndrome and 
assuming that columns indexed by /C U {jt} are erased. 
Compute e«,j, = e^f • 

c) Update the received array T and the syndrome array S 
by 

T{y,x) ^ T{y,x) - e^^jg-x^'y'^' 

Siy,x) <- S{y,x) - £«^,jj-Td_i(a;;QjJ-T„(?/;/3jJ . 

4) For every h £ (m), apply a decoder for Cgrs using row h 
of S as syndrome and assuming that columns indexed by /C 
are erased. Let E be the m x n matrix whose rows are the 
decoded error vectors for all h G (m). 

5) Compute the error array 

£^ ( H,j'Ea I H-'E^ | ... | U'l^E^-^ ) . 

Output: 

• Decoded array T — f of size mx n. 



Fig. 3. Decoding of errors and erasures of type (Tl), (T2), (T4), but not of 
type (T3). (See Section ITV-B] ) 

jo, Ji, ■■■,jw are all distinct, while = juj+i = •••== 
Furthermore, d and w should satisfy the inequalities 



w + T + p <d-2 . (34) 

(While the inequality in (1331 ) is already part of the requirements 
in Theorem |5] we need the second inequality so that ( fT3] l will 
hold. Specifically, the inequality in ( |34] | says that the number 
of erroneous columns does not exceed d—1. Observe that the 
inequality 2t + p < d — 2 and the inequality in ( [33l l, together 
imply the second inequality therein whenever m < d — p.) 

Without any loss of generality, we will also assume that 
^Ktjt 7^ for every £ e {■&). The set {je}ee{w+i) will be 
denoted hereafter by C. When i? = 0, we formally define w 
to be and C to be the empty set. 

Let the modified syndrome a be the unique m x 



matrix that satisfies 

a{y, x) = S{y,x) ■ ]^ (1 - a-jx) (mod x'^"^) , 

and let S be the m x {d—l—p) matrix formed by the columns 
of fj that are indexed by {p, d—1). We recall from (fTji that 
p = rank(^) = rank{{E)juc')- 

If p > 2w + 2, then we can regard the columns that are 
indexed by C as full block errors (namely, errors of type (Tl)), 
and the conditions of Lemma [5] will still be satisfied, namely, 
we will have 

2[t + w + 1) + p < d + p-2 . 

Therefore, we assume from now on that p < 2w + \. 

By ( fT3b we get that for every j e JlJC , column Ej belongs 
to colspan(5). In particular, this holds for j E C \ {jw}, in 
which case Ej (in polynomial notation) takes the form 

Ejiy) = EkJ • T,„(y;^«j) . 

Let the row vectors Uq, Oi, . . . , a,„_^_i form a basis of 
the dual space of colspan(S'), and for every i E {m — p), 
let ai{y) denote the polynomial of degree less than m with 
coefficient vector af, we can further assume that this basis 
is in echelon form, i.e., degao{y) < degai{y) < ... < 
degam_^_i(y) < to. This, in turn, implies that the degree 
of a{y) = gcd(ao(y),ai(y), . . . , a™„^_i(?/)) satisfies 

dega(?;) < p , 

namely, a{y) has at most p {< 2w + 1) distinct roots in F. 
Now, it is easy to see that for every ^ G F, the column 
vector (C'')fte{m) (also represented as Tm{y',£,)) belongs to 
colspan(S') (and, hence, to colspan(i?)j-u£'), if and only if 
^ is a root of a{y). In particular, /Jk^.j^- is a root of a{y) for 
every £ E (w) . We denote by TZ the root subset 

7^={(K,J) : a(/3,,,) = 0} , (35) 

and define the polynomial A(?;) by 

where rj = \TZ\. 

Consider the (rii-rj) x n matrix E = {eh,j)he{ni-r,),je{n) 
which is formed by the rows of A{y)E{y, x) that are indexed 
by {ipm). Specifically, 

n 

= X! ^h+-n-i,j , h e (m - 77) , j e (n) 

i=0 

(compare with (|32]|). Respectively, let S be the {rn—rj) x 
(d—l—p) matrix formed by the rows of A{y)S{y, x) that are 
indexed by {r],m). It readily follows that Ej^iy) = for 
£ e (w) and that 

(37) 
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Observe that the number of summands on the right-hand side 
of dJTl i is -d — w, and that number is bounded from above by 
771 — 2w — 1 < m — fj, < m — rj. This means that Ej^^ (y) ~ 
if and only if A(/3^^^^.^_ ) for all £ e {w, d). We also recall 
that 



rank(S') = rank{{E)ju 



(38) 



Next, we distinguish between the following three cases. 

Case 1: rj — fj.. By ( l38T l we must have Ej^ (y) = 0, which 
is equivalent to having A(/3~^^^^) = for all £ £ {•&). Thus, 
assuming this case, we have C CTZ, and the decoding problem 
then reduces to the one discussed in Section IIV-BI 

Case 2: 77 = /i - 1. If Ej^ (y) = then £ C 7^. Otherwise, 
it follows from (l38l l that each column in S must be a scalar 
multiple of Ej^. The entries of Ej^, in turn, form a sequence 
that satisfies the (shortest) linear recurrence 



i=0 



{K,j)eTZ' 



where 

7^ 



' = {{^e,J^) : i e (w,^) and A(/3-,i^.J ^ o} . 

Indeed, this recurrence is uniquely determined, since the num- 
ber of entries in Ej^ , which is m — 77 = m — /.j + 1 > m — 2w, 
is at least twice the degree \TZ'\ (< -d ~ w) of B{y). The 
recurrence can be computed efficiently from any nonzero 
column of S using Massey's algorithm II2TI . We now have 
£ C 7^ U 7^', where 

|7^u7^'| = \n\ + \n'\ 

< 7] + 1!} — W 

<2w + d-w 

<1} + W 

< m . 



So the decoding problem again reduces to that in Section lTV-BI 
Case 5; 77 < /i - 2. If Ej^ (y) = then (again) C C U. 
Otherwise, the conditions of Lemma |3] hold with respect to 
E and to the matrix Z formed by the rows of A{y)Z{y,x) 
indexed by {rj, m) (each such row is a codeword of Cgrs)- 
Hence, we can decode E. Next, we observe from (l37T i that for 
j = jw, the vector Ej{y) can be seen as a syndrome of the 
column vector 

Ke(m):A(/3-i)#0 

with respect to the following [m—rj) x 777 parity-check matrix 
of a GRS code: 



^GRS — 



where 



£{m—r]),K,G{m) ' 



/3l,A(/3-]) ifA(/3-^)=^0 
1 otherwise 



(39) 



(40) 



And since the Hamming weight of £* is at most — w < 
{m — rj)/ 2, we can decode £* uniquely from Ej (again, under 



the running assumption that j = j^)- Thus, for every k such 
that A(/3~^) 7^ 0, we can recover the error value and 
subtract it from the respective entry of T, thereby making 
TZ a superset of the remaining symbol errors. The problem 
though is that we do not know the index Therefore, we 
apply the above process to every nonzero column in E with 
index j ^ /C. A decoding failure means that j is certainly 
not jtu, and a decoding success for j 7^ j^j will just cause us 
to incorrectly change already corrupted-columns in T, without 
introducing new erroneous columns. We can then proceed with 
the decoding of T as in Section IIV-BI 

Figure |4] presents the implied decoding algorithm of a 
combination of errors of type (Tl), (T2), and (T3), provided 
that the type-(T3) errors satisfy the assumptions laid out at the 
beginning of this section; as said earlier, these assumptions 
hold when m < d — p and the number of type-(T3) errors is 
at most 3. Steps [lJ{3] |6al and |7]-[8] in Figure |4] are essentially 
applications of steps in Figures |2] and |3] Next, we analyze 
the complexity of the remaining steps in Figure |4] staring 
with Step Hi] A basis in echelon form of the left kernel of 
S can be found using 0{(Pm) operations in F, and from 
this basis we can compute a{y) using 777 applications of 
Euclid's algorithm, amounting to 0{dm^) operations. The 
set TZ can then be found in Step |4b] via Chien search, 
requiring O (77777 •min((i, 777)) operations, followed by 0{d^m) 
operations to compute the matrix S in Step|4c] The complexity 
of Step |5] is dictated by Step |5b] therein which, with Chien 
search, can be implemented using 0{m?n) operations. Finally, 
Step l6bl requires 0(^d{d + m)m) operations in F. In summary, 
the decoding complexity of the algorithm in Figure |4] amounts 
to 0[{d + 777)77777) operations for syndrome computation and 
the Chien search, and 0{d{d + m)m) operations for the other 
steps. 

Acknowledgment 
The authors thank Erik Ordentlich for helpful discussions. 

Appendix A 
Analysis for Example[T4] 

We show that (fTSl l implies the inequality (fTTI i. Specifically, 
we assume the following value of re6q{N,D) 

re6,{N,D) = riog,(V,(7V-l,Z?-2) + 1)1 , 

where 

D-1 



y,{N-l,D-2)^J2 



i=0 



N-1 



and show that, for 1< Di < D2 < ^p^N + 2, 

redq{N,D2) > red,(7V, A) + (i?2 - A) • 
To this end, it suffices to show that for 1 < D < + 1, 

red,(iV, D+1) > redg{N, D) + 1 . (41) 

For the above range of D we have 

N -i N-D + 1 q 



> 



> 



D-l ~ q-1 



l<i< D-l 
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Input: 

• Array T of size m x n over F. 

• Set IC of indexes of column erasures. 
Steps: 

1) Compute the m x (d— 1) syndrome array 



It follows that (for the above range of D), 



S = ( HoTo I HiTi 



Hn-lTn-l ) HJrS 



2) Compute the m x (d— 1— p) matrix 5* formed by the columns 
of a;) njgK(l ~ '^j^) 'hat are indexed by (p, d — 1). Let 
^ = rank(S'). 

3) (Attempt to connect assuming < /i/2.) Apply Steps [3]- 
|4]in Figure |2] (with = /C) to the modified syndrome array 
a{y, x), to produce an error array E. If decoding is successful, 
go to Step H] 

4) a) Compute the greatest common divisor a{y) of a basis of 

the left kernel of S. 

b) Compute the set TZ and the polynomial A(y) as in i35i - 
OU. Let ri = \n\. 

c) Compute the (m—rj) x (d— 1— p) matrix S formed by 
the rows of A{y)S{y,x) that are indexed by {r],m). 

5) If ri — fi — 1 then do: 

a) Compute the shortest linear recurrence B(i/) of any 
nonzero column in S. 

b) Compute the set 

7^' = {(a., j) : A(/3-j) / and B(/3-j) = O} . 

c) If |7^'| = degB(i/) and |7^'| < m - r? then update 

7^ ^ 7^ u 7^'. 

6) Else if 77 < /i — 2 then do: 

a) Apply Steps [2}E1 in Figure |2] (with A' = K.) to the 
syndrome array S, to produce an en'or array E. 

b) For every index j /C of a nonzero column of E do: 

i) Apply a decoder for the GRS code with the parity- 
check matrix Hq^^ as in l[39t-(l40t. with Ej as 
syndrome, to produce an eiTor vector £* . 

ii) If decoding in Step 6(b)i| is successful then let E* = 



Hj£j and update Tj ■<— Tj — £* and S{y,x) 
Siy,x)-E*(y)-Ta-iix;a,,). 

7) Apply Steps |2]^ in Figure [3] to 5, /C, and 7?,, to produce an 
error array E. 

8) Compute the error array 



£^ ( H^'Eo I ffr'£^i I •■• 

Output: 

• Decoded array T — £ of size m x n. 



Fig. 4. Decoding of en'ors and erasures of type (Tl), (T2), (T4), and 
with restrictions on errors of type (T3). (See Section HVCl ) For the sake 
of simplicity, we assume that there are no erasures of type (T4). 



which readily implies that 

^iV-l\ q (N-l 



> 



q-1 

and a simple induction on i thus yields 



1 < i < D-1 



^N-l\ ( Q V ( iV-1 



D-lJ - \q-lJ \D-l-i 



0<i< D-1 . (42) 



D-l 



,(iv-i,i^-2)^(,-i)--5:r:;_j(, 



< 



iV-1 
D- 



1 ('^-D'-^E'^-' 



1=1 



Namely, 



y,{N~l, + 1 < ('^_ J] {q 1)^-2 



Multiplying both sides hy q — 1 yields 



iV-1 



D-1 



(<7 - 1) • (V,(iV-l, D-2) + 1) < (^^_^ j {q - 1) 

and, so, 

<Z-(V,(iV-l,D-2) + l) 
< V,(iV-l, D-2) + 1 + (^^2 J) {q 1)^-1 

= Vg(iV-l,i?-l) + l . 

We conclude that 

1 + log, {y,{N-l,D-2) + l)< log, (V,(iV-l,i?-l) + l), 
and (HTt finally follows by taking ceilings on both sides. 

Appendix B 
Analysis for Example[T6] 



We derive the lower bound (125) on the expression 



^(r,-r,_i)(D,_i-l-2T) 



V 



E E 

1=1 j=ri_i+l 

2z? + 1 



^E 

j=i 

^E 

j=i 
1? 

= E 

i=i 



J 



2d + l 
3 

2i? + l 



"2i9 + l" 

2i?+ 1 
r^-i + 1 

2§+l 
U-i + 1 



- 1 



- 1 



- 1 



where the penultimate step follows from ru > 2^9. The lower 
bound ( I25] ) immediately follows by the known expression for 
harmonic sums ifTOl p. 264]. 
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